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ABSTRACT 

The Draco dwarf spheroidal galaxy (dSph), with its apparent immense mass to light 
ratio and compact size, holds many clues to the nature of the enigmatic dark matter. 
Here we present deep photometric studies of this dwarf galaxy, undertaken with the 
MegaCam Camera at the Canada-France-hawaii Telescope, the Wide Field Camera at 
the Isaac Newton Telescope and the Wide-Field and Planetary Camera on board the 
Hubble Space Telescope. The new photometric data cover the entirety of the galaxy, 
and reach i' = 24.5 at 50% completeness, significantly deeper than previous panoramic 
studies, allowing searches for tidal disturbances of much lower surface brightness than 
has been possible before. With these improved statistics, we find no evidence for 
asymmetric disturbances or tidal tails that possess more than 3% of the stars found 
within the centre of the galaxy. We find that the central stellar density, as probed by 
the HST data, rises into the central 0?5. Uncertainties in the position of the centroid of 
the galaxy do not allow us to determine whether the apparent flattening of the profile 
interior to 0f5 is reliable or not. Draco is therefore a flawless dwarf galaxy, featureless 
and apparently unaffected by Galactic tides. 
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1 INTRODUCTION 

The nature of the most common type of galaxy in the 
Local Group, the dwarf spheroidals, remains a mystery. 
Were these small galaxies formed early in the Universe, 
as predicted in the hierarchical galaxy formation model 
I White fc Reeslll978l ), and were many of them later incor- 
porated into the Milky Way? Though the ACDM model 
gives currently the best means to explain the formation of 
these galaxies, this scenario has the well-known problem 
of a discrepancy of 1-2 orders of magnitude between the 
number of low mass structures predicted by this theory and 
the observed number of satellite galaxies (|Moore et al.l 19991 : 



Based on observations obtained with MegaPrime/MegaCam, 
a joint project of CFHT and CEA/DAPNIA, at the Canada- 
France-Hawaii Telescope (CFHT) which is operated by the Na- 
tional Research Council (NRC) of Canada, the Institute National 
des Sciences de l'Univers of the Centre National de la Recherche 
Scientifique of France, and the University of Hawaii. 
Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained from the data archive at the Space Telescope 
Institute. STScI is operated by the association of Universities for 
Research in Astronomy, Inc. under the NASA contract NAS 5- 
26555. 



iKlypin et alj|l999h. Th eir apparent alignments in great cir- 
cles l|Lvnden-Belll 11976) have led s ome authors t o propose 
a tidal (non-cosmological) origin (Kroupa 1997), possibly 
by th e early interaction betw een Andromeda and the Milky 
Way (|Sawa fc Fuiimotoll2004l) . 

Thus the study of dwarf spheroidal galaxies is of impor- 
tance in our attempt to understand the formation of galax- 
ies in general and also the distribution of the dark matter. 
These small galaxies are most easily examined in the Local 
Group due to their proximity, and of the dwarf spheroidal 
galaxies in this region Draco is arguably the most interest- 
ing in terms of its importance in placing constraints on the 
distribution of dark matter on small scales. 

The Draco dSph is a small, faint, and metal-poor Galac- 
tic satellite situated about 71 kpc from the Sun (see Ta- 
ble 1). The stellar population is very old ( > 8-10 Gyr) and 
very metal poor ([Fe/H] = —1.8 ±0.2 dex) l|Apparicio et al.l 
1200 lib Its total i-band lum inosity is (L/Lq) = 2.4±0.5x 10 s 
l|Odenkirchen et al.ll200ll ), and King model (|Kindll962l ) fits 
to the stellar profile, show that it has a core radius of 7.7', 
and a tidal radius of 40.1'. Under the assumption of Virial 
equilibrium, the high stellar velocity dispersion implies an 
extremely high mass t o light ratio of about 146 ± 42 in So- 
lar units, according to lOdenkirchen et al.l d200ll) (he reafter 
O01) or 330 ± 125 according to iKlevna et all l|2002i ) (here- 
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Table 1. Properties of Draco 



Parameter Value Reference 

a (J2000) 17 h 20 m 13.2 s 001 
5 (J2000) +57°54'54" O01 



(l,b) (86. 4°, +34.7°) 

d 71±7kpc O01 

V r -293.8l2'7kms- 1 Harereaves et aJ. C 1996) 

cry km s — 1 llargreaves^t_aL (1996) 

R c 7.7' or 0.16kpc O01 

Rt 40.1' or 0.83 kpc O01 

E(B-V) 0.029 

[Fe/H] -1.8±0.2dex Apparicio et al. (2001) 

M 3.5 ±0.7 x 1O 7 M O01 

(M/L) ; 146 ±42 O01 

M v -8.8 Mateo (1998) 



after K02). Several other studies (e.g., iKlessen et al.ll2003l : 
lLokasl ! 2002). also conclude that it is probably a very strongly 
dark-matter dominated galaxy. 

In many of its properties Draco is very similar to the 
Ursa Minor dSph, another low-mass Milky Way satellite. 
Surveys of the latter galaxy, however, have shown the pres- 
ence of tidally-induced features wh ich attest to the non- 
equil i brium state of that galaxy llOlszewski fc Aaronsonl 
1 19851 ; IBellazzini et al.ll2002l ; IPalma et al.ll2003l ). and imply 
that the dark matter fraction may well be much lower in 
Ursa Minor than might be naively expected. In contrast, 
previous studies of Drac o have shown it to have a smooth 
symm etrical shape (O01: IPiatek et al ] |200ll ; IBellazzini et al] 
2002). Indeed, no study has yet found evidence for a dis- 
torted distribution of stars in the halo of the Draco dSph, 
which would be expected if it did not possess a massive 
dark matter mini-halo to protect the stellar component from 
the strong tidal field o f the Milky Way. Recently however, 
IWilkinson et al] (|2004l ) (hereafter W04) discovered a break 
in the light profile of Draco at ~ 25' and a sharp decline 
in the velocity dispersion at ~ 30'. They concluded that the 
outer part of Draco is filled with a dynamically colder stellar 
population, which could be caused by an external tidal field. 

Our aim in this work is therefore to search more thor- 
oughly for tidal perturbations in the halo of the Draco dSph, 
using the new generation of wide-field cameras. The outline 
of the paper is as follows: §2 describes the observations and 
the data reduction processes, §3 details the selection crite- 
ria for Draco stars, while §4 discusses the analysis of the 
photometric survey, the different star populations that are 
present, and the search for substructures in the outskirts of 
this small galaxy. Finally, §5 and §6 are dedicated to the 
discussion of our results and the conclusions of the study. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 MegaCam and INT surveys 

The Draco dSph galaxy was observed with the wide-field 
MegaCam camera at the Canada-France-Hawaii Telescope 
(CFHT) in service mode during the 2004A and 2005A ob- 
serving seasons. Megacam is a mosaic of 36 4612 x 2048 
EEV chips, offering a field of ~ 0.96 x 0.94 square degrees. 
We chose to observe seven fields, six around the galaxy and a 



comparison field four degrees in declination below the centre 
of the galaxy. Taken in photometric conditions, the images 
totalled 950 s, 1090 s and 1700 s per field in, respectively, the 
MegaCam g, r and i bands. The seeing ranged from 0"46 to 
1"12, with the majority of the images being taken in seeing 
better than 0"8. 

The images were pre-processed by the CFHT "Elixir" 
pipeline yielding debiassed, flat-fielded images, with i-band 
fringing artefacts removed. Photometric standards observed 
over the season are used by "Elixir" to measure the zero- 
point in each passband for every pointing. The resulting 
images were subsequently passed through the photometry 
pipeline developed for the Wide Fie ld Camera (WFC) o f 
the Isaac Newton Telescope (INT) (|lrwin fc Lewis! 1200 ll ). 
suitably altered to deal with MegaCam data. 

C atalogue gen eration follows the precepts outlined by 

llrwinl (| 19851 . Il996l ) and includes the facility to: automati- 
cally track any background variations on scales of typically 
20-30 arcsec; detect and deblend images or groups of im- 
ages; and parameterise the detected images to give various 
(soft-edged) aperture fluxes, position and shape measures. 
The generated catalogues start with an approximate World 
Coordinate System (WCS) defined by the known telescope 
and camera properties (the WCS distortion model) and are 
then progressively refined using the 2MASS (2 Micron All 
Sky Survey) astrometric catalogues to give an internal pre- 
cision generally better than 0.1 arcsec and a global external 
precision of 0.25 arcsec. These latter numbers are solely de- 
pendent on the accuracy of the astrometric catalogues used 
in the refinement. 

All catalogues for all CCDs for each pointing are then 
processed using the image shape parameters for morpholog- 
ical classification in the main categories: stellar, non-stellar, 
and noise-like. A sampled curve-of- growth for each detected 
object is derived from a series of aperture flux measures 
as a function of radius. The classification is then based on 
comparing the curve-of-growth of the flux for each detected 
object with the well-defined curve-of-growth for the general 
stellar locus. This latter is a direct measure of the integral 
of the point spread function (PSF) out to various radii and 
is independent of magnitude, if the data are properly lin- 
earised, and if saturated images are excluded. The average 
stellar locus on each detector is clearly defined and is used 
as the basis for a null hypothesis stellar test for use in classi- 
fication. The curve-of-growth for stellar images is also used 
to automatically estimate frame-based aperture corrections 
for conversion to total flux. 

Various quality control plots are generated by the 
pipeline and these were used to monitor characteristics such 
as the seeing, the average stellar image ellipticity (to mea- 
sure trailing), the sky brightness and sky noise, the size of 
aperture correction for use with the "optimal" aperture flux 
estimates (here "optimal" refers to the well-known property 
that soft-edged apertures of roughly the ave rage seeing ra - 
dius provide close to profile fit accuracy, e.g.. lNavlorlll998l ). 
The "optimal" catalogue fluxes for the g, r, i filters for each 
field were then combined to produce a single matched cata- 
logue for each pointing and the overlaps between pointings 
were used to cross-calibrate all the catalogues to a common 
system to within an accuracy ~ 1% across the survey region. 

For completeness we also present here an earlier 
panoramic dataset obtained with the INT WFC on the 
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nights of the 11th, and the 13th to 15th of October 2002. 
The WFC is a mosaic of four 4k x 2k EEV chips, offering 
a field of ~ 0.29 square degrees. A total of 11 pointings 
were observed in photometric conditions, with exposures of 
1800 sec in the equivalent of the Johnson V passband and 
1800 sec in the equivalent of the Gunn i band at each point- 
ing. The seeing ranged from 0"7 to l"l. The INT WFC im- 
ages were reduced in an identical manner to the MegaCam 
observations, with the exception that the pre-processing was 
performed with the INT pipeline. A small offset of 0.106 
magnitudes was applied to the i-band INT magnitudes to 
bring them into agreement with the MegaCam values. 

Tables and list the photometry resulting from the 
MegaCam and INT surveys, showing g, r, i-band MegaCam 
magnitudes (in the AB system), V and i-band INT magni- 
tudes (converted into the MegaCam AB system) , with their 
uncertainties, a classification parameter in each passband 
(c g , c r , a and cv), and the equatorial coordinates for each 
object. This classification index is a flag that indicates the 
likely nature of the objects. A classification index of 1 corre- 
sponds to a likely galaxy, — 1, . . . , —9 to a point-sources (in 
order of decreasing confidence of stellarity) , and to an arte- 
fact or noise. The INT photometry is only used to calibrate 
the HST V-band photometry. 



2.2 HST survey 

We also analysed a series of HST images of a field in the 
centre of the Draco dSph. The data were taken in 1995 and 
1999 with the Wide Field and Planetary Camera (WFPC2), 
in two spectral bands, V (F606W) and i (F814W). The total 
exposure time amounts to 4000 sec in F606W and 4500 sec 
in F814W. The wide-field (WF) chips of the mos aic were 
reduced with the DAOPH OT/ALLSTAR package (IStetsonl 
Il987l : IStetson et al.l [l990l ). A point spread function (PSF) 
prev iously constructe d for a study of the globular cluster 
M4 l|lbata et al.lll999h in each of the WF chips was used to 
fit the profiles of all detected sources in the WF mosaic. This 
yielded, for each detected object, the coordinates, the mag- 
nitude, a parameter describing the profile shape (sharpness) 
and the x value of the PSF fit. The equatorial coordinates 
of the sources were determined by comparison to the Mega- 
Cam survey. The instrumental WFPC2 F606W and F814W 
magnitudes were calibrated onto the INT V-band and Mega- 
Cam i-band systems, respectively, by assuming a simple off- 
set with no colour terms. The correspondence between the 
bands is reasonably acceptable, as we show in Fig. 1. 

The images listed in Table 5 do not cover exactly the 
same area of sky, and have offsets of up to ~ 80 pixels. Since 
we were primarily interested in measuring the stellar density 
profile in the inner regions of the Draco dSph, we retained 
only the overlapping region in common to all 14 exposures, 
thereby ensuring that the region was surveyed to a uniform 
depth. 



3 SAMPLE SELECTION 
3.1 Completeness 

The deep HST photometry in the centre of the galaxy is 
useful to help define the stellarity criteria and also to de- 



+ + + \ * 



t + 



19 



20 



21 



V (INT) 



22 



23 




(MegaCam) 

Figure 1. HST F606W photometry was calibrated as a simple 
offset from the INT V band (upper panel). Likewise, HST F814W 
photometry was calibrated with MegaCam i (lower panel). (The 
increased scatter at i=20 is due to photometric variability of hor- 
izontal branch stars). 



termine the completeness of the MegaCam survey. We find 
that adopting MegaCam sources with classification indices 
of either -1 or -2 in all three passbands results in zero false 
objects (i.e. zero sources not present in the HST survey). The 
completness for this selection is represented with a dashed 
line in Fig. 2. We also show, with a solid line, the complet- 
ness when we select all the MegaCam objects (classification 
index from 1 to -9). To i — 24 both selections have com- 
pleteness > 50%. Down to this limit, the adopted selection 
on classification indices does not make us lose a substantial 
fraction of stars. It should be noted that the incompleteness 
is affected to a significant extent by the crowding, which is 
maximal in the central region probed by the HST field, but 
will be much less problematic in the outskirts of the galaxy. 
Fig. 2 therefore represents the worst completeness over the 
survey. 

3.2 Star selection 

Fig. [3] shows the colour magnitude diagram (CMD) of the 
MegaCam stars in the central region of Draco within a radius 
of 6'. To reduce the contamination by foreground stars and 
background galaxies, we construct a selection region cover- 
ing the red giant branch (RGB) to sub-giant branch (SGB) 
on the central CMD. As the red giant branch has very low 
contrast compared to the "background" in some fields (or 
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Table 3. The first 10 rows of the INT catalogue. 
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Table 4. HST field IDs 



Figure 2. The completeness of the MegaCam dataset is com- 
pared to the HST survey for all the catalogue (solid line) and for 
the classification indices -1 and -2 (dashed line). 



name 


passband 


exposed time 
(s) 


date 


u2oc0101t 


F606W 


200 


1995 


u2oc0102t 


F606W 


1100 


1995 


u2oc0103t 


F606W 


1000 


1995 


u2oc0104t 


F814W 


200 


1995 


u2oc0105t 


F814W 


1100 


1995 


u2oc0106t 


F814W 


1100 


1995 


u5erl401r 


F606W 


400 


1999 


u5erl402r 


F606W 


400 


1999 


u5erl403r 


F606W 


500 


1999 


u5erl404r 


F606W 


500 


1999 


u5erl405r 


F814W 


400 


1999 


u5erl406r 


F814W 


500 


1999 


u5erl407r 


F814W 


500 


1999 


u5erl408r 


F814W 


500 


1999 



may simply not be present), this CMD selection will help 
reduce contaminants in the analysis below. A fairly broad 
selection is chosen here to ensure that we do not reject a 
significant fraction of members of the dwarf galaxy. 

Clearly a simple CMD filter is not an optimal way to 
select Draco-like stars, particularly towards fainter magni- 
tudes, where the discrimination between stars of the dwarf 
galaxy and background galaxies is difficult. We therefore try 
to implement a CMD filter that will enhance the signal of 



the dwarf galaxy by filtering out those CMD regions that 
degrade the total signal to noise. To this end we take the 
background field (Fig. 4) as representative of the contami- 
nation by foreground and background Galactic stars as well 
as point-like distant galaxies that have managed to evade 
the stellarity criterion. In general, the signal to noise ratio 
in a given CMD bin will be simply 



y/Sn% + pSn%' 
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Table 5. The first 10 rows of the HST catalogue. 
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25.304 


0.075 
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57 53 36.9 
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0.031 


0.760 


0.062 


24.806 


0.050 


0.950 


0.136 
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Figure 3. The colour-magnitude diagram of the central 6' of the 
Draco dSph, showing g magnitude as a function of g — i (left) and 
g — r (right). The initial CMD selection is indicated. 




Figure 4. As Fig. 3, but for the comparison field. 



where ur is the number of stars in the CMD bin in the re- 
gion under consideration, na is the number of stars in the 
comparison field, and / is the ratio of the geometrical ar- 
eas of the two fields. Assuming Poisson noise, this simplifies 
further to 

S/N= ^-f n * _ 

\/ n R + P n B 

Let us consider first the region inside 6' (displayed in Fig. 3), 
taking the stars within the selection polygon. We bin the 
CMD in intervals of 0.05 mag in g-i colour and 0.25 mag in 
g-band magnitude. For each bin, we calculate the signal-to- 
noise ratio, and we sort the resulting array. Starting with the 
CMD bin of highest signal to noise, we examine each new 
bin in decreasing order of S/N testing whether the additional 
bin adds to the total S/N, in which case the bin is accepted. 
It transpires that upon applying this process we accept the 
whole region within the polygon of Fig. 3, i.e. in the region 
within 6' of the galaxy centre, all of the CMD bins contribute 
usefully to the total S/N. 

However, this will not the case on the outskirts of the 
galaxy, where ur — friB ~ 0. In the outer regions the signal 
from the dwarf galaxy is of course extremely weak, so one 
cannot apply the same procedure as above. We use instead 
the CMD within 6' as a model. Now the signal to noise ratio 
becomes: 



\/riR + g 2 n D 

where n D is the number of Draco stars in the central (6') 
region and g is an unknown normalisation factor combin- 
ing the ratio of geometrical areas between the area of in- 
terest and the central region and the relative frequency of 
stars from the dwarf galaxy in the outer region compared to 
the central region. Since ur >> g 2 no in this background- 
dominated case, we find that S/N * Cn D (n R )-^ 2 , with 
C a constant of proportionality whose value we need not 
know. We proceed as above, binning the CMD, calculating 
the signal-to-noise ratio in each bin (to within the multi- 
plicative constant), and sorting the resulting array. We find 
that in the outskirts of the galaxy at radii beyond 45' an op- 
timal CMD selection includes the main-sequence turnoff and 
parts of the horizontal branch and blue-straggler regions. 

The left and middle panels of Fig 5 show the result 
of these two selection methods, where we present the full 
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Figure 5. The left and middle panels show, respectively, the 
result of the CMD filters constructed to optimise signal to noise 
in the central regions (SI) and in the outer regions of the galaxy 
(S2). The right hand panel shows the weight array of the matched 
filter analysis. 



dataset of point-sources, filtered by the two CMD filters. 
These two CMD filters allow us to construct a star selection 
appropriate for the centre of the dwarf galaxy (SI, on the 
left-hand panel), and the other (S2, in the middle panel) 
that enhances the signal to noise in the external parts of the 
galaxy, where the background is high. We retain 25670 stars 
for the SI selection and 5385 stars for S2. Fig. 6 shows the 

spatial distribution of the two selections. 

The matched-filter technique (see e.g., iRockosi et al.l 

I2002T I provides an alternative means to enhance the signal 
of an expected CMD structure in the presence of contami- 
nation. We apply this procedure also to our survey, taking 
the stars within the polygon of Fig. 3 to be representative of 
the desired signal, with a background sampled by the com- 
parison field shown in Fig. 4. The matched-filter method 
assigns a relative weight to each star, according to its CMD 
location (the map of the matched-filter weights is displayed 
in the right-hand panel of Fig. 5). In this way it is more 
efficient than the methods used above to select the samples 
"SI" and "S2" since one can make use of the entire dataset. 
However, the drawback is that one must bin all quantities 
of interest, which is not always convenient. We therefore use 
all three methods, applying each where appropriate. 

Fig. [7] shows the selection boxes used to investigate the 
spatial distribution of stars at other evolutionary phases 
which are clearly present from visual inspection of the CMD: 
the blue and red horizontal branch (hereafter BHB and RHB 
respectively), and the blue stragglers (hereafter BS). These 
boxes are defined from the (g-i,g) CMD. The spatial dis- 
tribution of these selections are shown on Fig. [8] where we 
display maps for the BHB, RHB and BS populations, which 
contain 605, 2542 and 973 stars, respectively. 



4 ANALYSIS 

4.1 Stellar density profile 

We next construct the stellar density profiles for the two 
selections. However, there are two practical problems that 




£ (arcmin) 

Figure 6. The distribution of sources selected from the opti- 
mized CMD selection regions displayed in Fig. \E\ the external 
optimisation S2 is shown on the upper panel, whereas the central 
optimisation SI is portrayed on the lower panel. The positions of 
the sources are shown in standard coordinates (gnomonic projec- 
tion). 
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Figure 7. The selection boxes for the BHB, RHB and BS evolu- 
tionary phases, superimposed on the photometry of stars drawn 
from within the central 6'. 
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Figure 8. The spatial distributions of the three selections, BHB, 
RHB and BS, arc shown. 



need to be solved. First, the survey is riddled with holes 
caused by CCD artefacts around bright stars, so these unus- 
able areas of the fields need to be masked out. We eliminated 
circular regions around the positions of bright USNO and 
Hipparcos stars which have magnitudes B < 16, choosing 
appropriate magnitude-dependent radii to cover the holes 
in the MegaCam survey. The second problem is the gaps 
between CCDs. The chipsets of the MegaCam camera are 
distributed on four rows, with two large gaps. These gaps 
(which are clearly visible in Fig. 6) are sufficiently large to 
influence measurements of the stellar density profile. 

To trace the stellar density profiles, we construct a suc- 
cession of ellipses with different major axis values (through- 
out this work, we use s to denote an elliptical radius, and 
r to denote a circular radius). We choose the best shape 
parameters found by O01, a ellipticity equal to 0.30 and a 
position angle of the major axis of 88°. For each optimised 
selection, we construct the surface brightness profile, taking 
into account the holes made by the bright stars and the gaps 
between CCDs, and remove a background level. The back- 
ground count level is determined from the comparison field. 
In Fig. 9 we show the resulting profiles before background 
subtraction (upper panel) and after background subtraction 
(lower panel). 

We next fit structural models to the SI and S2 profiles 
and compare our results with O01 and lWilkinson et al.ll20o3 
(herea fter W02). We use three different models, the King 
model (|Kindl 19621 ): 



So 



1 



y/i + (S/Rc) 2 y/l + (Rt/Rc) 2 J 

the exponential model (jSersidl 19681 ): 

E(s) = E exp(-( S / J R ) n ) 

and the Plummer model: 

E(s) = E i?o/(l + (s/Roff . 

The fitting is performed by \ 2 minimisation, with the 
ellipticity and the position angle of the major axis taken as 




Figure 9. The stellar profile of the Draco dSph galaxy. The upper 
panel shows the profiles prior to background subtraction, while 
in the lower panel the mean count level from the comparison field 
has been removed. The black histogram shows the profile of the SI 
sample, the S2 sample is displayed in red, and the profile resulting 
from the matched-filter analysis is shown in green. 



constant. The results are summarised in table 6, and the 
corresponding profiles are displayed in Fig. [TO] 

The extensive spatial coverage of our survey provides us 
with information far from the core of Draco, up to ~ 100'. 
In contrast to W04, who found a break at 25' in the light 
profile with the INT dataset, the profiles on Fig. [9] based on 
either the S2 selection or from the matched-filter analysis do 
not show evidence for a break. We suspect that the presence 
of the break in the W04 study could have been an artefact 
of background subtraction, which is very difficult at these 
low count levels. 

The fits of Fig. 10 give a visual impression that the data 
can be represented quite well by these simple models. How- 
ever, the reduced x 2 values listed in Table 6, are very poor, 
and even the best of the models (the exponential model fit 
to the Matched-Filter sample) can be rejected with better 
than 99.99% confidence. Thus the data are now of such ex- 
cellent quality that these simple models no longer provide 
satisfactory fits. We suspect that dynamically self-consistent 
triaxial models will be required to reproduce the structure 
of this dwarf galaxy, however that is beyond the scope of the 
present article. 

We also trace the stellar density profile for the BHB, 
RHB and BS selections. Fig. [TT] shows the three resulting 
profiles. A fit with a Plummer model yields scale radii of 
ll'.8 ± 3!8, 7!6 ± 0!3 and 8!l ± 0!2 respectively for the BHB, 
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Table 6. Model fits for the parameters of the King, exponential, 
and Plummer profiles using the SI, S2 and Matched-Filter (MF) 
samples. The quantity X 2 l(N — f) gives the reduced-x 2 for the 
fit. We also indicate the values found in earlier studies. 



King 


data 


rc(') 


rt(') 


X 2 /(N~f) 




SI 
S2 
MF 
001 


8.15 ±0.02 
7.18 ±0.05 
7.63 ± 0.04 
7.7 ±0.2 


45.0 ±0.4 
48.9 ±0.9 

45.1 ± 0.6 
40.1 ±0.9 


5.35 
4.61 
3.31 


Exponential 


data 


ro(') 


n 


x 2 




SI 
S2 
MF 
O01 


8.07 ±0.05 
6.80 ±0.11 
7.43 ±0.10 
7.6 ±0.1 


1.20 ±0.01 
1.07 ±0.01 
1.15 ±0.01 
1.2 ±0.1 


3.97 
3.42 
2.27 


Plummer 


data 


ro(') 




x 2 




SI 
S2 
MF 
W02 


10.73 ±0.01 
10.09 ±0.02 
10.06 ±0.01 
9.71 ±0.1 




9.55 
5.03 
8.01 
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Figure 10. Fits of a King (red), an exponential (green), and a 
Plummer (blue) model to the stellar density profiles derived from 
the Matched-Filter sample (upper panel) the S2 sample (middle 
panel) and the SI sample (lower panel). The error bars represent 
lcr uncertainties. 




Figure 11. The density profiles of the BHB (blue dashed line), 
RHB (red dotted line) and BS (black solid line) selections. 



RHB and BS selections. We find that widening the BHB 
box by allowing redder horizontal-branch stars into the se- 
lection results in a profile that is closer to that of the RHB. 
Thus the paucity of BHB stars does not allow us to confirm 
the visual impression from Fig. 8 that the BHB stars are 
more dispersed than the R HB and BS stars. C onsequently, 
Draco, unlike the Scu lptor (jTolstov et alj|2004l ) and Fornax 
|Battaglia et aljfeood ) dSph, does not possess a BHB pop- 
ulation that is significantly more extended than the RHB 
population. Whether this attests to a difference in early 
star-formation or simply to the low statistics of BHB stars 
is currently not clear. 



4.2 Smoothness of the stellar density profile 

The stellar profiles displayed in Fig. 10 appear ap- 
proximately smooth, but are they really featureless 
apart from their monotonic decrease? To investigate this 
point we implement a varia nt of a test devised by 
iKuhn. Kocevski fc Fleckl (|2004l ) (hereafter KKF), who in a 
preliminary study of the stellar profile of the Draco dSph, 
found the very interesting result that the inner regions 
within 20' are not as smooth as one would expect from Pois- 
son noise statistics. Following KKF, we divide our Matched- 
Filter sample into 1000 radial bins, using the same ellipti- 
cal coordinate as in Fig. 10. Each bin contains 24 stars to 
maintain uniform statistical significance. However, since our 
survey suffers from gaps between CCDs and holes due to the 
haloes of bright stars, we calculate the density pi in each bin 
i, and also the density uncertainty 8 pi, taking into account 
the Poisson noise in each bin, and the Poisson noise in the 
background subtraction. 

In a similar way to KKF, we take groups of 5 bins at 
radius R t , and calculate the reduced \ 2 of the group, test- 
ing the null-hypothesis that the density is constant over the 
small radial interval. To this end, we define: 



3=0 



T>) 2 



where p is the weighted mean density in the group of 5 
bins. Large values of x\ would indicate large deviations from 
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Figure 12. The colour-magnitude distribution of stars in the 
WFPC2 field. A CMD selection polygon is applied to these data 
to remove a small number of probable contaminants. The inset 
shows the spatial location of these stars at the centre of the dwarf 
galaxy. 




Figure 13. The stellar profile of the HST sample (uppermost 
points), is compared to the profiled derived from the matched- 
filter MegaCam sample, previously displayed in Fig. 9. The deeper 
HST data give better statistics in the galaxy centre, allowing us 
to probe the density profile down to s = 0'.3. 



Poisson noise, and if many such bins were found, it would 
imply significant small spatial-scale variations in the den- 
sity. However, contrary to KKF, we find that in the region 
within R < 20' the observed distribution of \\ does not 
differ significantly from the expected distribution; indeed 
a Kolmogorov- Smirnoff test indicates that the probability 
that the observed distribution is drawn from the expected 
X 2 distribution is 2%. The hypothesis that the stellar profile 
of Draco is smooth on small-scales therefore cannot be con- 
fidently rejected. Since any defects in the survey that were 
not accounted for in our processing will have the effect of 
driving down this probability, one should interpret this test 
as evidence that the stellar profile is smooth. 

4.3 The central profile of Draco 

We include the HST survey in this discussion of the spatial 
properties of the Draco dSph, as the resulting deep pho- 
tometry complements the wide-field MegaCam data, giving 
much better statistics in the centre of the galaxy. The CMD 
derived from the point-sources in the much deeper WFPC2 
data is shown in Fig. 12. (We define point sources to be 
those objects in Table 5 that have \ < 1-5 anu sharpness 
between and 0.2 in both the V and i-bands). The main 
sequence turnoff and part of the RGB and HB are clearly 
visible. In this central region of the dwarf galaxy, the number 
of contaminants to the point-source catalogue is very small, 
but we choose to implement a simple CMD selection box (as 
shown in Fig. 12) around the stellar populations of Draco 
to avoid the obvious contamination. The selection inside the 
polygon contains a total of 1122 stars. 

We now trace the surface density profile of this selec- 
tion of HST stars. We divide the sample into 11 bins each 
containing 102 stars, and as with the MegaCam dataset, 
we compute the stellar density taking into account the gaps 
between the CCDs. Fig. 13 compares the stellar density pro- 
files of Draco obtained with the HST data to the MegaCam 
matched-filter selection. In Fig. 14 we show more clearly the 
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Figure 14. The stellar density profile of the central regions of 
the Draco dSph, as measured from the HST survey. Each bin has 
102 stars. For comparison, with the dashed line we display the ex- 
pected projected density profile of an NFW model of Virial Mass 
10 s Mq and concentration c = 28.6 (the density normalisation 
has been chosen to minimise the x 2 fit to the data). The dotted 
line shows the profile of an NFW model with c = 4.5. The full line 
shows the exponential model fit previously in the upper panel of 
Fig. 10 (but suitably normalised to minimise x 2 )- 



region probed by the HST data; the higher source density re- 
sults simply from the greater photometric depth. The mean 
star-count level in the HST survey is 296 ± 41arcmin~ 2 , 
however the probability that these density are constant with 
radius is 0.8%. Thus the profile derived from the HST star- 
counts requires that the stellar density in the inner regions 
of the galaxy continue to increase towards the centre. 

In cosmological simulations, dark matter halos follow 
a "univer sal" profile, like th e NFW (|Navarro et al.l Il997r ) 
or Moore l)Moore et al.ll 19991 ) models. Though these models 
work well on large scales such as those of galaxy clusters, 
the question about the form of the inner regions of low-mass 



10 Mathieu Segall, Rodrigo Ibata, Michael Irwin, N., F. Martin & Scott Chapman 



dark matter halos is not re solved: are their centres cuspy or 
cored (|de Blok et aUl200ll )? Though the distribution of the 
luminous matter need not follow that of the underlying dark 
halo, it is nevertheless interesting to compare the observed 
distribution of stars to the expected dark matter models. 

In Fig. 14 we compare the observed central stellar num- 
ber density profile to a cusped profile model, that we choose 
to parametrise with a p rojected NFW model of mass 10 8 M© 
l)Lokas fc Mamonll200ll ). We adopt a normalisation to min- 
imise x 2 ! s ° that the only free parameter is the concentra- 
tion c of the NFW model. The mean concent ration of a halo 
of this mass is c = 28.6 (M accio et alJ [20061 ): this model is 
shown with a dashed line, and has a probability of only 0.1% 
of conforming to the data. However, there is a wide range 
in the concentration o f cosmological halo s, with a typical 
scatter of 0.4 in log(c) (|Macci6 et alj|2006h . The dotted line 
shows a fit for c = 4.5, an extreme (2a) value for halos of 
mass 10 8 Mq . The probability that the observed profile can 
be drawn from such a model is 8%. 

The marginal discrepancy between the NFW models 
and the data stems solely from the innermost point of the 
survey. We have checked this datum carefully, and could find 
no artificial reason why it should be substantially lower than 
the trend. This may reveal a real departure due to a change 
in the mass to light ratio of the central regions, though it 
may also simply point to the possibility that the centre of 
the dwarf galaxy is not precisely at the location measured 
by O01. Their centroid resulted from a fit to the data over 
the entire galaxy, so it is reasonable that there could be 
an offset in the very central 0!3. Indeed, this is perfectly 
consistent with the quoted uncertainty of 0!3 in the O01 
centroid position. The MegaCam survey, due to the large 
gaps between CCDs, is not well suited to improve on the 
O01 centroid measurement, and neither is the HST survey 
presented here due to its very limited spatial coverage. If 
we omit the innermost datum, the probability that the HST 
profile is drawn from a c = 28.6 or a c = 4.5 model is 48% 
and 14%, respectively. 

Nevertheless, it is interesting in this context that the 
simple exponential model that we fit earlier to the Matched- 
Filter sample provides a good representation to the data, 
including the innermost datum. This is shown with the full 
line in Fig. 14, where we have simply scaled the previous fit 
to best match the observed profile. The probability that the 
data are drawn from this exponential model is 26%, while 
rejecting the innermost point yields a probability of 28%. 

4.4 Search for tidal features 

One of the particularities of the Draco dSph based on pre- 
vious studies of this system is its undisturbed appearance. 
To try to place better constraints on the absence of tidally- 
induced features we count the number of stars at radii be- 
tween 45' < R < 120' probing beyond the King model tidal 
radius to the end of the survey. Each annulus is divided 
into 36 equal parts, and we count the number of stars in 
each annulus segment. In Fig. 15 we show the background- 
subtracted star-counts and their la uncertainties in each an- 
nulus segment for the S2 selection. None of these 10°-wide 
bins exceeds la, so we conclude that no significant overden- 
sities are detected. 

To test the sensitivity of this method, we make fake tidal 
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Figure 15. Counts in azimuthal bins around the Draco dSph 
(the error bars represent la uncertainties) for stars beyond 45'. 
The data are drawn from the "S2" selection, which is appropriate 
for the outer regions of the survey. 



tails by redistributing a fraction of the stars from the central 
region of the galaxy along the direction of the major axis of 
the galaxy. For every randomly drawn star from within a 
20' radius of the galaxy centre, we randomly reassign the 
coordinate along the major axis of the galaxy so as to have 
a constant density stream over the length of the surveyed 
region. With a tidal tail representing 3% of the selection 
of the central stars, we are able to detect the tail at 4.9<r 
(summing two bins around 0°) and 3.3a (from summing 
two bins around 180°). This result indicates that a tidal 
tail consisting of only 3% of the central body of Draco is 
detectable by this method. 

The absence of any other deformities in the outskirts of 
the galaxy is further corroborated by the stellar density map 
in Fig. 16, which has been constructed from the matched- 
filter sample. The uncertainty contours at the 1, 3 and 5<r 
levels are also shown. These data clearly show that Draco is 
a very regular and symmetric system. 



5 DISCUSSION 

5.1 The inner density profile 

The HST survey presented above shows clearly that the in- 
ner stellar density profile increases inwards until at least 
0!5. Inside that radius, the profile may flatten or decline, 
or we may simply have adopted a slightly incorrect cen- 
troid. The current data are not good enough to discrim- 
inate between these possibilities. CDM simulations show 
dark matter profiles with a power-law slo pe between -1 and 
-1.5, as shown by numerous recent studie s (| Colin et al ] |2004l ; 
iGhigna et ai1l2000l ; lHavashi et alJbOOSt ). However, low sur- 
face bri ghtness galaxies app ear to have flat central mass 
profiles l|de Blok et al.ll200lf ), in stark contradiction to the 
theoretical predictions. These authors also show that the 
profiles are better fit with a model with a power-law index 
of -0.2 ±0.2. 

A flat central density profile has also been deduced for 
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Figure 16. Matched-filter map of the MegaCam survey region, 
shown on a logarithmic scale and over the same spatial region as 
in Fig. 6. The lines show significance contours of the signal to 
noise ratio at 1, 3 and 5o\ 



the Ursa Minor dSph, from an a nalysis of the surviv ability 
of a morphological substructure l|Klevna et al.l I2003T) . They 
found that the stars of this substructure are dynamically 
colder than the stellar population of the dwarf galaxy, con- 
cluding that the substructure was a stellar cluster accreted 
by Ursa Minor, and whose stars now orbit unbound inside 
the halo of the dwarf galaxy. They demonstrated that this 
type of scenario is only possible if the dark matter profile 
around Ursa Minor has a flat core, as otherwise the stars 

rapidly lose their coherence. 

Some numerical simulations (jRicotti fc Wilkinson! 

|2004 > can produce dark matter halos with flat cores in 
agreement with those observations. In their cosmological 
simulations, iRicotti fc Wilkinson! (|2004T ) created this type of 
halo at z > 10, by taking into account the effect of radiative 
pressure. They are able to produce galactic systems similar 
to the Local Group. To make this, they selected at z — 10 
galaxies with cored dark matter halos. Then they evolved 
without any peroid of fusion. The main result is the core of 
each dwarf galaxy remains flat and stable until z = 0. The 
modelled velocity dispersion profile was foun d to have the 
same shape as those found observationally bv lKlevna et alj 
(2002). Another method to obtain dark matter halos with 
flat core was proposed by iRead fc Gilmord (|2005l ). They 
discovered it is possible to flatten the core of a potential 
by two short periods of mass- lo ss, separated by a phase 
of reaccretion of gas. Recently, iMashchenko et all (|2006l ) 
showed that the CDM cusp may be removed for halos in 
the 10 7 to 10 10 M range by including stellar feedback into 
the simulations. 

Clearly our starcounts are not directly measuring the 
dark matter profile; the stellar profile will have resulted 
from a complex interplay between baryons and the dark 
matter. With complementary kinematic measurements we 
could hope to use the Jeans equations to disentangle the 
two. However, though samples of s everal hundred s t ars ar e 
now being obtained in this galaxy IWilkinson et al.l l|2004h , 
the HST photometry we have analysed above probes the in- 
teresting very inner parts of the system, where kinematics 



are as yet unavailable in sufficient numbers to allow a mea- 
sure of the velocity dispersion, let alone its gradient. (The 
paucity of kinematic data is due to the low surface number 
density of bright giant stars needed for spectroscopic obser- 
vations). Nevertheless, a clear flattening of the starcounts 
in the central regions of the galaxy away from an NFW 
model would have lent support to the presence of a core in 
the Draco dSph. Similarly, if the stars were found to follow 
an NFW profile into the very central regions it would call 
into question the possibility that the dark matter is cored 
in this dwarf galaxy. Unfortunately, the uncertainty in the 
position of the centroid of the galaxy does not allow us to 
draw a clear conclusion on this matter. If we do not include 
the innermost datum, the HST data follow an NFW model 
closely. Including the innermost datum forces us to accept 
a model with a flat inner profile like the exponential model 
shown in Fig. 14. Further observational work is required to 
measure accurately the position of the centre of the dwarf 
galaxy, and to improve on the statistics of the starcounts in 
the central few arcmins. 

5.2 Absence of tidal effects 

We have searched for tidal disturbances in the outer parts of 
the Draco dwarf spheroidal galaxy, using a stellar selection 
designed to optimise the signal to noise ratio in regions of low 
contrast. However, no evidence was found for the presence of 
an extra-tidal population, and we limit the fraction of such 
a component to less than 3% within the spatial extent of the 
MegaCam survey. This result, in conjunction with the find- 
ing that the surface brightness profile is consistent with that 
of an unperturbed dynamical system, suggests that Draco 
has led a quiet existence. This finding also lends support 
to models in which the stars are shielded by a dark matter 
mini-halo from the tidal forces of the Milky Way. 

A similar conclusion about the non-detection of tida l 
streams in Draco has been made by iKlessen et all i|2003h . 
who studied the properties of blue horizontal branch stars 
in Draco with the SDSS dataset. They choose this stel- 
lar population because it is distributed all over Draco and 
the contamination by other astronomical objetcs is mini- 
mal. Their surface brightness profile of this stellar popu- 
lation does not possess any break, in agreement with our 
much deeper MegaCam results presented here. They com- 
pleted this study with a numerical simulation of the dwarf 
galaxy, without a dark matter component and affected by 
the tidal effects of the Milky Way. By comparing the mor- 
phology of the simulation to the observations of the hori- 
zontal branch of Draco, they concluded that such a purely 
baryonic model is not compatible with reality. Their numer- 
ical simulations show a wide horizontal bran ch that is not 
observed for Draco. Thus lKlessen et all i|2003h conclude that 
Draco is not a renmant of a tidally disrupted galaxy, but a 
dwarf galaxy with a substantial dark matter halo. 

The fact that we find no tidal stream is not alto- 
gether unexpected. iMaver et all (|2002h conducted cosmolog- 
ical simulations to attempt to model the behaviour of dwarf 
galaxies with tidal debris, such as the Carina dSph. They 
found that to reproduce this type of galaxy, the dark mat- 
ter halo must have a flat core or else it can possess a steep 
inner profile but with very low concentration (< 7). If the 
concentration is higher, the orbits must be highly tangential 



12 Mathieu Segall, Rodrigo Ibata, Michael Irwin, N., F. Martin & Scott Chapman 



to reproduce the observations. In the case of dwarf galaxies 
with high mass-to- light ratio, like Draco or Ursa Minor, they 
showed that it is difficult to create large tidal streams. Only 
small tidal features could be formed, but at low contrast. 



6 CONCLUSIONS 

We have presented a study of the structure of the Draco 
dwarf spheroidal galaxy, using photometric data obtained 
with the MegaCam camera on the CFHT telescope and the 
WFPC2 camera aboard the Hubble Space Telescope. These 
two datasets give us a new detailed view of this small galaxy. 

Our deep photometry allows us to probe further than 
previous studies into the heart and the outskirts of the dwarf 
galaxy. We find that the central stellar density increases in- 
wards from ~ 2' to 0'.5. We find reasonable agreement if 
we choose to model the central HST stellar profile with a 
cusped model, parametrised as a projected NFW model of 
mass 1O 8 M0. Inside 0!5 the profile appears to flatten out, 
or drop, though this may be an artefact of adopting an in- 
correct value for the position of the centre of the galaxy. 
However, we find that an exponential model gives an equally 
good description of the data. 

The derived stellar profile is very smooth on large scales 
and does not appear to have a break, contrary to some previ- 
ous claims. We also do not find evidence to support a claim 
that the stellar distribution is clumpy on small scales. We fit 
the stellar profile with various simple models, and improve 
slightly the model parameters compared to earlier studies. 
However, neither a King model, an exponential or a Plum- 
mer model give satisfactory fits to the observed profile over 
the entirety of the galaxy. The data depart slightly from each 
of these idealised models, and we are forced to reject them 
all at high confidence. However, this discrepancy is more a 
testament to the outstanding quality of the data. A more 
refined model will have to be developed to properly model 
these data; this is however beyond the scope of the present 
paper. 

We also search for possible tidal substructures in the 
outskirts of the dwarf galaxy. No evidence for the presence 
of a tidal stream is found. Monte-Carlo tests show that if any 
such structure is present it cannot represent more than 3% of 
the main body of the system. Thus we find that Draco is an 
undisturbed system, most likely in dynamical equilibrium. 
The Draco dSph therefore continues to be an excellent target 
for studies of the dark matter on small scales. 
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